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PREFACE 


The method of calculating the electromagnetic pulse (EMP) produced 
by a high-altitude burst, as presented in this Note, was originally 
discussed at the Joint EMP Technical Meeting-~First Annual Nuclear EMP 
Meeting, held at the Air Force Weapons Laboratory in September 1973, in 
a paper by the author entitled A Nwnerical Example of the Effect of 
Atmospheric Seattering on Predicted EMP Environments, Several requests 
for a more detailed version of the method were received from attendees 
at the meeting. During the following years a number of detailed com— 
parisons were made by the author and the Air Force Weapons Laboratory 
(AFWL/DYT} of calculated electromagnetic pulse environments from speci- 
fied bursts using the author's method and the AFWL's CHEMP code. When 
the same input parameters were incorporated into the calculations, the 
numerical agreement in the pulse amplitude versus time-was excellent-- 
typically within about 25 percent at times of principal interest, i.e., 
to times well past the time of maximum pulse amplitude. 

This Note should be of interest to analysts involved with EMP envi- 
ronmental predictions and also to those desiring additional physical in- 
sight into the generation of the EMP signal and how it is influenced by 
the gamma source time characteristics, by atmospheric scattering, and by 
residual ionization of the atmosphere. It is being published at this 
time to make available to those concerned with EMP environments an al- 
ternative and physically descriptive method for determining EMP levels. 
It contains the unclassified general theory and analysis portions of a 
more extensive classified report of the same title published in 1974, 

The Note was prepared under the sponsorship of the Defense Advanced 
Research Projects Agency's study project "Future Strategic Communica- 


tions Systems and Nuclear Survivability Issues." 
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SUMMARY 


This Note presents the essential details of a three-dimensional 
method of calculating the radiated electromagnetic signal that is 
caused by magnetic deflections of a distribution of Compton electrons 
produced in the upper atmosphere by prompt gamma radiation from a 
high-altitude nuclear detonation or a series of detonations. The 
method differs from that previously developed for such calculations 
in that the solution is obtained from summing the radiation fields 
from the individual electrons in a three-dimensional volume instead of 
combining the individual electron motions to determine a time— and 
space-dependent current from which the radiation field is evaluated 
using a one-dimensional approximation to the solution of Maxwell 
equations. 

In the Note the effects of the time characteristics of the source 
gamma (y) output, the effects of atmospheric scattering in reducing 
the coherent radiation from the Compton electrons, and the effect of 
residual atmospheric ionization (preionization) in reducing the obser- 
ver signal amplitude are developed and illustrated numerically. 

The three-dimensional characteristics of the source are brought 
out in illustrative numerical examples. These examples give additional 
physical insight into the source characteristics and demonstrate the 
importance of scattering and atmospheric ionization, when present, in 


influencing the EMP amplitude/time characteristics. 
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I. INTRODUCTION 


The purpose of this Note is to outline and numerically illustrate 
the essential details of a fully three-dimensional method of calcu- 
lating the radiated electromagnetic signal from magnetic deflections 
of a distribution of Compton electrons produced in the upper atmos- 
phere by prompt gamma radiation from a high-altitude nuclear detona- 
tion. The three-dimensional method evaluates the radiation fields (at 
a given time) from individual electrons. It uses Lienard-Wiechert 
potentials and sums the fields from all the spatially distributed elec- 
trons which contribute to the observer field at the given time. Rep= 
etition of the process for various times after the detonation leads 
to the full determination of the signal amplitude as a function of time. 
The method outlined--development of the solution in terms of the 
radiation field from an individual electron-~is different from the 
method which has been generally employed for calculation of the radia- 
tion. _In previous solutions, “1? the method has been to combine the 
motions of all Compton electrons to determine a time- and space- 
dependent electric current and then to calculate the radiation from 
this current. Both methods, properly applied, must give the same 
answer for the radiation field from a given Compton electron configura- 
tion, since both are derived from Maxwell's equations. A formal proof 
(2) 


of this has been given by Sollfrey for the zero conductivity case. 
As will be shown, summing the fields of the individual electrons 
requires relatively little computational effort and readily permits 
the inclusion of atmospheric scattering and absorption effects in the 
solution. Atmospheric scattering was not considered in Ref. 1 but has 
been shown °3°4) to be, in general, an important factor in determining 
the radiation field. Also, the present method provides additional 
physical insight into the process of generation of the radiated elec- 
tromagnetic pulse (EMP) arising from the magnetic deflection of 
Compton electrons. 
Section II of this Note gives details of the principal steps of the 


electron field summation method of EMP calculation. For simplicity in 


presenting the method, vertical geometry and a magnetic field perpen- 
dicular to the burst/observer path are used, and numerical values are 


illustrated for all steps in the development of the method. 


II. DETAILS OF THE SUMMED ELECTRON-FIELD METHOD 


A. GENERAL DISCUSSION 

Figure 1, assuming vertical geometry, shows the geometrical 
aspects of the summed electron-field method. A short burst of y radia- 
tion is emitted by a source well above the altitude where the emitted 
y radiation effectively interacts with the atmosphere. The y radiation 
impinges on the atmosphere below the burst and produces Compton elec- 
trons which deflect in the earth's magnetic field and radiate an 
electromagnetic signal which travels to the observation point. The 
electromagnetic pulse is determined at an illustrative location directly 


below the burst by appropriately summing the time-dependent electric 


fields produced at the observation point by all the electrons in the 
ellipsoidal volume that can contribute at a given time. In the 
calculation, the time-dependent effects of atmospheric scattering on the 
field produced by the electron are explicitly included, as is the 
reduction in the field caused by absorption by secondary ionization, 
which exists between a given electron and the observation point. 

The first signal reaches the observer at the time R/C (C = velocity 
of light) after the first y radiation from the burst. This signal 
results from Compton electrons produced on the direct path between the 
burst and the observer. The ellipse, much shortened in Fig. 1, illus- 
trates the boundary of the ellipsoid of revolution about the burst/ 
observer axis that contains all the radiating electrons that can 
contribute to the observer field at a time t - R/C, or t, after arrival 
of the first signal. It is easily shown that the length of the axis, 

X, of the ellipsoid at an axial distance ry from the burst and ty from 


the observer is 


x =(——=* (ep) 
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Fig. 1 Geometrical factors for 100 km burst and ground observer (vertical path) 


To illustrate the dimensions of the problem, let 


R=r,+r, = 10° meters 
t= 107° seconds or 1 shake 


=7*x 10* meters 


ial 

tr, = 3x 104 meters 
Then 

X = 355 meters 


or any time t, in shakes, 
X = 355 Vr meters. 


Consider a disc whose radius is X, as shown in Fig. 1. Clearly, at 
times up to about 107° seconds (or 100 shakes), the dimension X is very 
small relative to the distance of the disc from the burst or to the 
observer. Hence the y flux from the burst can be considered constant 
across the disc. ; 

The y flux at a given height and the number of first scattered 
Compton electrons produced at this height are readily calculated by use 
of total mass attenuation coefficients for air, such as those given by 


(5) (6) 


Evans, and a standard reference atmosphere. For any arbitrary 

y spectral distribution, the y spectrum can be divided into several 
small segments with an appropriate average y energy per segment. The 
Compton electron production for each segment is evaluated by means of 
the procedure of the moncenergetic calculation, The total Compton 
production is then the sum of the Compton electrons produced by the 
several segments. Figure 2 illustrates the height distribution of 
Compton electrons produced by a prompt y source of 0.01 KT yield. The 
figure portrays vertical geometry, a burst height of 100 km, and a 
source emitting 1.6-MEV monoenergetic y radiation. Also shown on Fig. 2 
are the height distributions of total secondary electron density 
production (N,) and the electron density production due only to the 
first generation Compton electrons a). 


The Compton electrons produced by Compton scattering of y radia- 


tion of a given energy have initial energies and directions as given 


Compton electrons /meter/km 
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Fig. 2 — Compton electron production vs height 
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by the Klein-Nishina formatas:* Figure 3 illustrates the electron 
energy and relative number per unit angle as a function of angle from 
the burst/observer axis for 1.6-MEV y radiation. For observer times 
{i.e., true signal time at the observer location) of principal interest 
(a few shakes) and typical EMP geometry such as that in Fig. 1, the 
initial distributions shown in Fig. 3 are essentially constant through- 
out the ellipsoid, since the angle between the line from the edge of 
the disc to the observer and the burst/observer axis is small. In the 
analysis in this Note the initial distributions shown in Fig. 3 are 
assumed to apply throughout the Compton electron source region, confined 
to the ellipsoidal region of Fig. 1. This assumption results in an 
overestimate of the radiated signals that is negligible at early times 
but increases with time. 

Rather than deal with the radiation from electrons that have number 
versus angle and energy distributions as shown in Fig. 3, it is conve- 
nient in: both the summed electron-field method and the continuum-current 
method of Ref. 1 to replace the distributions by equivalent forward- 
directed electrons that will produce the same initial field intenstty 
at the observation point as is produced by the distribution of electrons. 

Since the initial number and the energy distributions of electrons 
are known from the Klein-Nishina formulas, one may precisely evaluate 
the normalized initial field produced by this distribution at the 
observation point by summing the fields of all the electrons in the 
distribution. The field, at a distance r, due to a relativistic 
electron with angle 6 to the direction of the observer and with a 
velocity 8C perpendicular to a magnetic field B is given by Lienard- 
Wiechert potential theory (see any text on electrodynamic theory; 


e.g., Ref. 7) and in simple form is 


2 2,4 
e“8 (cos 6 - 8) B (1 - 8) (2) 


E= 3 
4nme (Cr (1 - 8 cos 9) 


where e is the electron charge, ™ is the rest mass, ey is the per- 


mittivity of free space, and C is the velocity of light. 
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Fig. 3—-Energy and number distribution for Compton scattering: 1.6-MEV y 


The field is normal to the line between the electron and the 
observer in the direction of the electron’s acceleration caused by the 
perpendicular component of the steady magnetic field. Figure 4 shows 
the distribution of the contributions to the observed radiation field 
as a function of the initial angle 6 of the Compton electrons (as in 
Fig. 3). For Fig. 4, a magnetic field perpendicular to the direction 
of the y ray was used, and the fields were evaluated for all electrons 
in the cone at an angle 8; the area under the curve is proportional to 
the total field contributed by the distribution. Again using the 
Lienard-Wiechert relation, we can readily calculate the velocity or 
energy of that number of forward-directed electrons which will produce 
exactly the same initial field intensity at the observation point as 
the ensemble of Klein-Nishina distributed electrons. For 1.6-MEV y 
radiation, the equivalent forward electron is determined by this 
process to have a 8 of 0.87 or an energy of 0.53 MEV. This value is 
to be contrasted with the average Compton electron energy of 0.8-MEV 
which, if used, greatly overestimates the observer field. 

Figure 5 shows a plot of the energy of equivalent forward-directed 
electrons versus incident y energy for a range of y energies. Calcu- 
lations for Fig. 5 were performed as illustrated above for 1.6-MEV 
incident y radiation. 

The above analysis applies to the field intensity of electrons 
before they have been rotated by the magnetic field and gives precisely 


(8) 


the correct result for electrons of zero age. Sollfrey has extended 
the above analysis by comparing as a function of time the observer 

field intensity due to the rotated ensemble of Compton electrons with 
that from the rotated equivalent forward-directed electrons. His 
analysis shows that, for the times of primary interest in EMP, the 

field versus time due to the rotated equivalent electrons is practically 
the same as that due to the rotated ensemble of Compton electrons. Thus 
the approach of using equivalent forward-directed electrons is quite 
accurate and provides a significant reduction in the computational 
effort. Again it should be emphasized that the energy of the equivalent 
forward electron for producing the correct observer field is signif~ 


icantly smaller than the average electron energy of the Compton electrons. 


Relative E field due to electrons at angle & times 
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Fig. 4 — Relative radiated field vs angle for Klein=Nishina 
distributed Compton electrons 
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This is true for either the summed electron-field approach or the 
continuum-current approach in EMP caiculations. 

Having determined (1) the height distribution of the first- 
scattered Compton electrons which results from a given y impulse and 
specified burst/observer geometry and (2) the equivalent energy of 
forward-directed electrons that produce the same field as the Klein- 
Nishina distribution of Compton electrons, we can proceed to determine 
the observer field arising from the full three-dimensional distribution 
of Compton electrons. 

Figure 6 shows, for a thin layer of Compton electrons in the 
height region where they are produced by the y radiation from the 
source (as in Fig. 2), the burst/observer axis and the outer boundary 
of the disc from which the field of radiating electrons can contribute 
to the observer radiation field at a time 1, where t is simply the 
difference between travel time at the velocity of light from the burst 
to the observer along the burst/observer axis and travel time along the 
path to the outer edge of the disc. As was discussed above, the 
Compton electrons, which have a distribution of energies and directions, 
can be replaced for purposes of determining the radiation field at the ~ 
observation point by equivalent initially forward-directed electrons 
of appropriate energy. For the moment, to aid in understanding Fig. 6, 
consider the y radiation to be emitted instantaneously, i.e., as a delta 
function. At the edge of the disc, the electrons that contribute to 
the observer field at time tT have a velocity parallel to the burst/ 
observer axis (since the dise radius, X, is very much smaller than the 
burst/observer distance, shown on Fig. 1). Electrons along the burst/ 
observer axis have been turned because of acceleration by the magnetic 
field and have acquired a direction determined by the amount of turning 
in a time t. The electrons between the axis and the outer edge of the 
disc have been turned for times intermediate between 0 and 1 and have 
velocity directions as illustrated. Consider the intermediate path 


wots 


to electron "q" shown by the dashed lines in Fig. 6. The relative y 


radiation travel time to "q" plus the relative travel time of the 
signal radiated by the electron to the observer is less than t, say 
wan 


t - At. Thus the electron "q’' has rotated for a time At when it makes 


its contribution to the observer field at time t. 
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Fig. 6 — Sketch illustrating electron directions which 
produce observer field at time t 
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It is straightforward to calculate all needed information about 
the initially forward-directed electron of known energy in a specified 
magnetic field as a function of time after its birth. The needed 
information is simply the variable parameters involved in the Lienard- 
Wiechert potential relation. Since the electrons have initial velocities 
near the velocity of light, the calculations must be made relativisti- 
cally, by properly relating electron time and motion to observer time. 
Consider again electron "q.'' To determine its contribution to the 
observer field at observer time t, we need the Lienard-Wiechert param- 
eters for the electron at time At, where 4t is in observer time. In 
general, we need the Lienard-Wiechert potential for observer times 
from zero to whatever length of time the EMP signal characteristics 
are to be calculated. Such calculations are easily done by computer. 
Figure 7 shows the normalized Lienard-Wiechert potential of an electron 
of velocity 0.87C (appropriate for a 1.6-MEV y source) moving in a 
Field of 0.6 gauss (perpendicular to the electron velocity vector) as 
a function of observer time. Note that, after 1.5 shakes, the field 
contributed by the electron is negative compared with the field from 


Q to 1.5 shakes. 


B. ILLUSTRATIVE CALCULATION USING A DELTA-FUNCTION Y SOURCE AND 
IDEAL CIRCULAR-ORBIT ELECTRONS 

Using Fig. 7, we can easily determine the observer field due to 
the thin layer of Compton electrons illustrated in Fig. 6. Consider 
Fig. 8, in which the abscissa is the distance from the burst /observer 
axis and the ordinate is the normalized contributions from electrons 
in the disc to the observer field. Again, to illustrate the approach, 
the electrons are considered to have ideal circular orbits (atmospheric 
scatter, absorption, and energy conservation effects on the ideal 
solution will be considered subsequently) and a delta-function source 
of y radiation is assumed. Figure 8 illustrates the situation for two 
observer times arbitrarily taken as 1 shake and 4 shakes for a burst 
altitude of 100 km and for the Compton electrons in a thin layer at 30 
km altitude. To the observer, the total field contributed per unit 


density of Compton electrons in the disc is simply 


Relative observer E field per electron 


= he 


1.6-MEV ¥ photon 
8B =.87 
B =.6 gauss 


i 2 3 4 5 6 


Electron rotation in shakes of observer time 


Fig. 7 — Lienard-Wiechert relative E field 
vs electron rotation time 
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Ep * 29X5E, f ae ax (3) 
() oo 

or 2nX, (maximum field intensity per electron) (area under the RE/X E, 
curve). The maximum field intensity per electron, Ey is simply the 
synchrotron field as determined by the Lienard-Wiechert expression 
(Eq. (2)) for the electron directed toward the observer (6 = 0). For 
8 = 0.87, B= 0.6 gauss, and x = 30 km, E = 4.28 X 10-* v/m/electron. 
For the 4-shake disc in Fig. 8, the area under the XE/X JE, curve is 
51.2 unit-meters. Thus the observer field at 4 shakes contributed by 
the layer of Compton electrons at 30 km is 27 X 710 m X 4.28 x 10714 
v/m/electron x 51.2 m= 9.8 xX 10? 


of Compton electrons). From Fig. 2, for example, at 30 km the Compton 


volt-meters/electron yx (density 


electron density produced by a 0.01-KT source of y radiation is shown 
to be 1.3 xX rot electrons/m“/km. Thus, the contribution to the 


observer field from the thin layer at 30 km is 


9 


9.8 xX 10 ” volt-meters/electron x 1.3 x 1921 electrons /m2/km 


E/h 


or 


E/h = 1274 v/m/km, where h is the layer thickness in kn. 
Following the same procedure for the 1l-shake disc of Fig. 8, we obtain 
1294 v/m/km. 

By repeating the. above procedure at other altitudes we can deter- 
mine the contribution per unit height at the desired time throughout 
the region containing Compton electrons. Finally, the total field at 


the observation point at a given time is 
H 


E= f[ E/h dh (4) 
) 


where H is the burst height. 

Figure 9 shows the results at a time of 4 shakes for the geometry 
of Fig. 1 and a y yield of 0.01 KT. The field intensity as calculated 
by the procedure leading to Fig. 9 is the ideal upper limit field 


Voits/meter/ilometer 
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Fig. 9 — Observer signal contributions versus height at 4 shakes 
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producible at 4 shakes by the source, geometry, and magnetic field 
used. If atmospheric scattering and absorption due to ambient atmos- 
pheric ionization were negligible, for the low y output used, the 
numerical value for field intensity would be reasonably close at the 
time shown (4 shakes) to the correct field, since other neglected 
factors are of minor consequence; however, addition of atmospheric 
scattering effects alone, as will be shown, can greatly reduce the 
field calculated for ideal circular orbits. Thus, practical meaning 
should not be attributed to the calculated field of Fig. 9. Meaningful 
values considering all necessary inputs to the calculation are presented 
later. The above material, it is believed, has adequately described 
the essential mechanics of the summed electron-field approach. In 
examination of Fig. 8, we see, for example, that at 4 shakes the 
observer field arises primarily from electrons in an annulus near the 
outer edge of the 4-shake disc. Actually, in the portion of the disc 
between the burst/observer axis and about 0.8 of the dise radius, : 
radiation from the electrons tends to reduce the observer field. This 
physical insight into the nature of the EMP source is a valuable 


feature of the summed electron-field approach. 


C. TIME-DEPENDENT Y SOURCE CALCULATIONS 
We next examine the effect on the above solution for ideal 

circular-orbit electrons of a source that emits y radiation over a 
period of many shakes. For convenience in illustrating the effect of 
a time-dependent source, we will assume an output rate that increases 
linearly to a maximum at 2 shakes and decreases linearly to zero 
between 2 and 8 shakes. We will outline the approach for handling a 
time-dependent y source and illustrate the approach by calculating the 
observer field at 4 shakes using such a source. Consider again Fig. 6. 
Since the time for which the observer field is to be evaluated is less 
than the time during which y radiation is being emitted by the source, 
at all distances from the burst/observer axis, electrons will be found 
having deflections appropriate for electron ages between zero and 4 
" shakes or less. At the outer edge of the disc the electrons have zero 


age but their number is also zero. At some intermediate position on 
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the disc (such as "q" in Fig. 6), the electron will have ages between 
zero and At (as is explained on p. 12). Thus, to obtain the field 
contributed by electrons at, say, "q", one must convolve the y rate 
output-versus-time curve with Fig. 7 to obtain the effective field 
contributed by the electrons of various ages. Figure 10 illustrates 
this process for locations on the disc where the excess travel times 

at the velocity of light between the source and the observer are 1 and 

3 shakes. At these locations the maximum ages of electrons contributing 
to the observer signal at 4 shakes are 3 shakes and 1 shake, respec- 
tively. The effective E source, at any location in the disc, is readily 
determined by the procedure illustrated in Fig. 10. At any location 

on the disc in terms of the abscissa time of Fig. 10, it is simply the 
ratio of the area under the convolved curves (labeled a xX b) to the 
total area under the relative y rate curve (b). Thus at 3 shakes in 
Fig. 10 the ratio is 0.138. Figure 11 shows for an observer signal 

time of 4 shakes the relative contribution to E as a function of disc 
delay time as obtained by this procedure. Clearly the procedure is 
applicable to any arbitrary y output rate-versus-time relation. 
Graphical convolution can be quickly and accurately performed. 

Figure 12 shows, for the same conditions as Fig. 8, the parameters 
determining the radiation from the thin disc at 30 km for the time- 
dependent y source illustrated in Fig. 10. Using the calculations 
previously given in connection with Fig. 8 for an observer time of 4 
shakes, we can readily determine the disc contribution for the time- 
dependent y source illustrated. The field is merely the field calcu- 


lated in Fig. 8 multiplied by the ratio of the area under the XE/X DE, 
33.1 
52.1 
of Fig. 8 with Fig. 12 brings out the difference in the physical 


curves of Fig. 12 and Fig. 8. This ratio is 


or 0.635. Comparison 


electromagnetic signal generation regions at the same observer time for 
a delta-function y source and for the time-dependent y source. In 
Fig. 12 the principal contributing region is about 500 m,from the 
burst/observer axis, and for Fig. 8 the principal contributing region 
is at the edge of the disc. The other differences are apparent. 

The total field intensity at the observation point is determined 


by repeating the procedure for Fig. 12 at several heights in the source 


SBxPYs | PUD | 4O BdsNOS anJoaJJa UIDJGO OF 
PJaly UOsjDaja YsIM ayns yndyno A Burajoauos yo uolosysnpyf—oQy{ “Sty 


@D4N0s yo aypa yndyno-4 satspjas (q) 


SHIGIO JOjNI41D ‘uCM2a[9 /3 a@AyD[es (p)——+, 


pup ani 4 aaippyay , 


i=) a foo] ~~ oO ” 7 ie) N - 
q . . : ° . > . . 


Woujd9j@/pjaly 9 aAlyojes 


Relative E 


-22- 


corresponds to edge of 4=shake disc 


Time = 0 
= 4 corresponds to burst/observer axis 


Time 
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Fig. 11 — Effective electron radiation vs time for 
time-dependent y source of Fig.10 
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region defining a curve similar to that of Fig. 9 for the delta-function 
source. Results of this procedure are not detailed here; however, the 
v/m/km curve versus height is essentially the same as that given in 

Fig. 9, except that the amplitude is uniformly reduced by a factor of 
0.635. Thus the total observer field at 4 shakes is essentially 0.635 

of that calculated (Fig. 9) for the delta-function source, or 17,000 v/m 
for a total y output of .01 KT. Again, practical meaning should not be 
given this calculated field value, since omitted factors, subsequently 


considered, greatly reduce the actual expected signal. 


D. EFFECTS OF INCLUSION OF ATMOSPHERIC SCATTERING OF COMPTON ELECTRONS 
Next we will proceed to a more meaningful solution by including 
the effect of atmospheric scattering in the field calculations. In- 
stead of the individual electrons following ideal circular orbits and 
producing a high degree of coherence in the summed field contributions, 
the electron directions rapidly become gaussianly distributed about the 
circular track they would follow. if not scattered. The mean squared 
angle of scatter, 8, is proportional to the distance, D, traveled by 


* 
the Compton electron after birth and is given approximately by 


92 ~0:5D 6 (5) 


8 20 
(Eypy) ~° 
where D is in meters 
o is the density at the electron altitude 


9, is sea level density. 


If the electrons were indeed initially forward-directed with an 
equivalent energy, Evry? as proposed on p. 8, scattering effects could 
be much more easily evaluated than in Compton scattering, where the 
electrons have a distribution of energies and initial directions rela- 


tive to the direction of the incident y radiation. Figure 13 illustrates 


* 

Fermi?) uses a factor of 0.7 instead of 0.5, but for electron 
energies near 1 MEV, the original work of Ref. 10 indicates the value 
0.5 as more reasonable. 
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in a general way the importance of scattering in deflecting the electron 
from the ideal circular orbit and, hence, in reducing coherence. For 
Fig. 13, electrons are considered to be initially forward-directed 
with an energy of 0.8 MEV (8 = .92). At 20 km, for example, the rms 
scatter angie becomes large relative to the magnetic turning angle, 
Pay at very early times. For the typical parameters illustrated, the 
electron, in the absence of scatter, contributes a positive field to 
the observer for about 1 shake, whereas with scattering effects 
included, the rms scatter angle reaches the zero field value in 0.1 
shake at 20 km. This does not mean that the field is zero in 0.1 shake, 
since the fields of all the electrons in the gaussian distribution will 
still contribute a positive value; however, it indicates the nature of 
the effect. We have approximated the full.scattering effect for various 
altitudes and electron parameters using both the continuum=current 
approach and the summed electron field approach. 

By summing the electron fields of all the electrons in the gaussian 
distribution using the Lienard-Wiechert potential theory discussed on 
p. 8 we have derived the results of Fig. 14 for 1.6-MEV y radiation 
and a perpendicular magnetic field of 0.6 gauss.’ Figure 14 shows the 
relative E field contributed by the average Compton electron (in the 
Klein-Nishina distribution) as a function of time for altitudes of 
20, 30, and 40 km. Also shown is the relative E field per electron 
based on the parameters of the circular-orbit forward-directed electron 
discussed on p. 10. Figure 14 illustrates the severe effect at 20 km 
of scattering in reducing the effective electron radiation and also 
shows that the scattering at 40 km is relatively small but not negli- 
gible. It is clear that neglect of scatter in EMP calculations can 
lead to sizable overestimates of the calculated fields, as originally 
shown in Refs. 3 and 4. 

Results derived for the relative E field, as a function of electron 
age for the various altitudes of interest, can be incorporated into 
the analysis directly. If the y source were a delta function, instead 
of using Fig. 7 for obtaining the dise distribution of Fig. 8, we 
would use the scatter curve for the appropriate altitude as illustrated 


in Fig. 14. For a time-dependent source, it is necessary to convolve 
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Circular orbits at all altitudes 


Relative electron E field 
2 
on 


0 — 
O 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.2 1.3 1.4 
Retarded time (shakes) 


Fig. 14—Relative electron field versus time with scattering 
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the relative E curve (with scatter) for the altitude of interest, such 
as those shown in Fig. 14, with the y rate curve (as was illustrated 
in Fig. 10 for the case of ideal circular orbits) in order to determine 
the contribution to the observer field from a given disc location at a 
given time. Figure 15 shows the results of convolving the relative E 
field of the ensemble of scattered Compton electrons with the y output 
rate illustrated in Fig. 10. It shows (as a function of time after 
the first y radiation arrives at any point in the y deposition region) 
the relative E field contributed to the radiation field by the average 
scattered electron in this region. At tT = or, the electrons contribute 
their maximum radiation field (or a relative field of 1.0) at all 
altitudes, with or without scatter. New Compton electrons continue to 
appear for a duration equal to that of the source pulse, and electrons 
of increasing age contribute decreasing radiation depending on the 
altitude of the disc. Thus the net result is a curve which decreases 
from a value of 1.0 to zero over a time somewhat longer than the dura— 
tion of y radiation from the source. Figure 15 shows the average 
Compton electron contributions that have been made at the times shown; 
at 2 shakes, for example, for the pulse illustrated in Fig. 10, only 
0.25 of the Compton electrons have been produced, and these electrons 
at 20 km produce an average field per electron 0.086 of the initial 
field from the forward~directed electron. 

Using results as illustrated in Fig. 15 and the procedure detailed 
for Fig. 8, we can proceed to evaluate the upper limit of the observer 
field considering a time-dependent y source and including atmospheric 
scatter. Still not included in the calculation at this stage are 
absorption effects and energy conservation. Figure 16 shows the dis- 
tribution of the E field source from the center to the edge of the l- 
and 4-shake thin discs (i.e., the discs at 30 km that contribute to 
the 1- and 4~shake observer E field) previously illustrated in Fig. 8 
for a delta-function source and with scattering neglected. Table 1 
details the calculated parameters used in plotting Fig. 16. In com- 
paring Figs. 8 and 16, note that the vertical scale of Fig. 16 is 0.1 
the scale of Fig. 8. Also plotted in Fig. 16 is the result for no 
scatter, previously shown in Fig. 12 for the 30-km thin layer contribution 
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Circular orbits 


Relative E field per electron 


0 2 4 6 8 10 
Time after pulse front (shakes) 


Fig. 15—Convolved scatter (Fig. 14) and time-dependent ¥ source (Fig. 10) 
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to the observer field at 4 shakes. The area under the XE/X EB, curve 
is proportional to the intensity of the observer field in units of 
v/m/km. The ratios among areas in the three cases considered in 


Fig. 8, Fig. 12, and Fig. 16 are as follows (at 4 shakes): 


Fig. 16 - Scatter and time-variable y source = 0.41 
Fig. 12 - No scatter and time-variable y source ’ 


Fig. 16 - Scatter and time-variable y source = 0.26 


Fig. 8 - No scatter and delta-function y source 


Table 1 


NUMERICAL VALUES FOR THE DISC AT 30 km CONTRIBUTING TO THE RADIATED 
E FIELD AT THE OBSERVER LOCATION AT 4 SHAKES OBSERVER TIME 
(t = 0 is instant of arrival of first signal) 


a b c a e £ g 
Radial 
Distance 
Path Delay from 
to Time after Burst/ Fraction 
Position Arrival of Observer Relative of Compton 
on Disc Gamma : Front Axis Convolved E Electrons E/E XE/K_E 
(shakes) (shakes) (meters} (per electron) Made fo) ome) 
0 4.0 0 +066 +672 +0443 0 
0.5 3.5 251 -083 578 .0483 O171 
1.0 3.0 355 -108 479 0517 .0259 
1.5 2.5 435 137 -370 +0507 +0311 
2.0 2.0 502 - 180 +250 -0450 -0318 
2.5 1.5 561 -246 141 +0338 -0267 
3.0 1.0 615 +335 +0625 +0209 -O281 
3.5 5 664 -52 +0156 -0081 -0076 
4.0 0 710 1.0 i¢] 0 ie) 
NOTES: 


Column f = Colum e X Column d 
Column g = Colum f X Colum c x 1/X, 


Xs 709 meters (Eq. (1)) 


=92= 


Using the same numerical factors previously derived for the 30-km 
layer (p. 17), we obtain 1274 X 0.26 or 331 v/m/km for the 30-km disc's 
contribution to the observer EMP signal at 4 shakes in observer time. 

By repetition of the above procedure, the 30-km disc contribution 
at 1 shake is determined to be only 19 v/m/km compared to 1294 v/m/km 
for the delta source and neglecting scatter (Fig. 8). In the time 
interval 0 to 2 shakes, it is apparent that the contribution of the 
disc to the observer field increases closely proportionally to (time)? 
With a delta-function y source, the signal increases with (time)!, and 
with a linearly increasing time-variable y source the number of 
contributing electrons increases with (time)?; thus the signal amplitude 
is proportional to (time)?. At 1 shake after arrival of the first 
y photons for the y time output rate assumed, as shown in Table 1, only 
0.0625 of the first-scatter Compton electrons which are ultimately 
produced (i.e., in 8 shakes) have been produced on the burst/observer 
axis and none have been produced at the edge of the contributing 1-shake 
disc. The above time characteristics are the same throughout the 
height range of the radiating Compton electrons; thus the total EMP 
signal produced by integration of the contribution of all the thin 
layers throughout the height region containing Compton electrons has 
the same amplitude-versus-time behavior at early times. For a given 
total y yield from the source, the rise characteristics of the 
radiated EMP depend, in general,on both (1) the rate of increase of 
the Y output rate and (2) the total time over which a given y output 


is produced. 


E. EFFECTS OF INCLUSION OF ABSORPTION EFFECTS FROM RESIDUAL 
IONIZATION AND SELF-IONIZATION 


General 

Ionization between the radiating electron sources and the observer 
will reduce the radiation, as was previously calculated in this Note 
and illustrated in Fig. 16. In the summed electron-field method of 
calculating EMP signals, the effect of this intervening ionization is 
calculated by use of standard propagation theory and analysis procedures. 


The field contributed by a given electron is reduced by the integrated 
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attenuation in the path between the electron and the observer at the 
proper time. Two sources of ionization, must, in general, be included: 
(1) ionization which existed in the region between the electron and 

the observer prior to the creation of the Compton electrons whose 
radiation is to be determined and (2) secondary ionization produced by 
the Compton electrons by inelastic collisions with the atmospheric 
gases. The first source, preionization, can be considered independent 
of time for times of interest in EMP problems. The second source, 
self-ionization, has an involved time-and-space dependence, which must 
be evaluated. 

Figure 17 shows an expanded view of Fig. 1 for the region between 
30 km and the observer. Boundaries of the ellipsoids for differential 
path delays between the source and the observer of from 0 to 4 shakes 
are shown. 

Consider again the radiation from the thin layer at 30 km which 
contributes to the total observer field at 4 shakes, as illustrated in 
Fig. 16 for the case of a time-dependent y output and including atmos- 
pheric scatter effects on the Compton electron radiation. Electron 
radiation from the edge of the 4-shake disc that contributes to the 
4-shake observer field traverses path "a" in Fig. 17 from the 30-km 
height to the observer; hence the reduction in signal contributed at 
the observation point from the edge of the disc is determined by the 
integrated differential absorption along path "a." Preionization will 
be essentially constant for EMP times of interest at a given altitude 
in any practical situation. Self-ionization at a given altitude is 
determined by the length of time Compton electrons have existed at the 
altitude. Along path "b", the burst/observer axis, Compton electrons 
have existed for 4 shakes independent of height below 30 km. Along 
path "a", however, the maximum time Compton electrons have existed and, 
hence, have been producing secondary electrons, is much less. For 
example, path "a" penetrates the 25-km height at a point where the path 
delay (from burst to observer at the velocity of light) is about 3.1 
shakes. Thus the signal contributed to the observer field at 4 shakes 
from the edge of the dise encounters self-ionization accumulated over 


4 - 3.1 (or 0.9) shakes, but along the burst/observer axis, the radiation 
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from electrons contributing to the observer field at 4 shakes encounters 
self-ionization, which has built up for 4 shakes. As a further illus- 
tration, consider path "c." Radiation from the 30-km disc contributing 
to the observer field at 4 shakes is attenuated by the self-ionization 
created in 2 shakes at the 30-km height and about 4 - 1.6 (or 2.4) 
shakes at 25 km. The differences in absorption at a given time for 
radiation from various portions of the disc may appear to cause compu- 
tational complexity; however, it will be shown to be relatively easily 


handled. 


Absorption Due to Self-Ionization 
To determine absorption produced by self-ionization along any path 


from the disc to the observer we first determine the electron density 
as a function of altitude and time after the arrival of the y radiation 
front. At the given altitude and time, ionization at all portions of 
the disc is the same, since the disc radius is much smaller than the 
distance from the burst to the disc; the time (in the observer's frame 
of reference) at which ionization builds up to a given value increases 
with radial distance on the disc from the burst/observer axis, as 
illustrated by the path delay curves in Fig. 17. Since the time 
history of first-scattered Compton electrons produced at a given alti- 
tude is known from a prescribed weapon y yield and y time history, the 
magnitude of the secondary ionization produced by the Compton electrons 
approximated can be as a function of time, considering the effects of 


(11) G2) i.e., the nonzero 


_ atmospheric scattering and ionization lag, 
time required for Compton electrons to produce low-energy secondary 
electrons. In Ref. 11, the problem was simplified by assuming that the 
high-energy electrons were initially foward-directed and that the 
electron energy was consumed in producing secondary ionization; to 

date no secondary ionization-versus-time calculations appear to have 
been published which use the full Klein-Nishina distributions of 
electron energy and electron direction as a point of departure and in 
addition account for energy loss to the electromagnetic field system. 
Such calculations, even for monoenergetic y radiation, are considerably 


more complex numerically than calculations using "equivalent" forward- 
'P y & q' 
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directed electrons. Figure 18, based on material presented in Ref. 11, 
shows the fraction of secondary electrons produced as a function of 
observer time for a forward-directed 0.742-MEV electron at 20-km 
altitude. At times of principal interest, i.e., at less than about 

8 shakes (since scattering has essentially randomized the Compton 
electron source for periods this long for altitudes below 40 km, as 
shown in Fig. 15), Fig. 18 may be applied to other altitudes between 
about 15 and 35 km by scaling by a factor fogs where 9 is relative 
density at altitude h, by including a correction for ionization lag as 
given in Ref. 12, and by scaling with energy as given by Evans. (9) 

To obtain the fractional secondary-electron density at a given 
altitude for a time-dependent source of y radiation, we must convolve 
the fractional ionization versus time (as illustrated in Fig. 18 for 
20 km altitude) with the source y rate-versus-time curve (as illustrated 
in Fig. 10). Figure 19 illustrates the convolution of Fig. 10 and 
Fig. 18 and shows the result for the fraction of secondary electrons 
produced at 20 km versus time for the y output rate of Fig. 10. The 
fraction applies to the total of secondary electrons produced by the 
first-scatter Compton electrons, which are shown as the curve Ny in 
Fig. 2. 

By the procedure illustrated in Fig. 19, we can determine the 
density of secondary electrons produced by Compton electrons as a 
function of time for any altitude in the y deposition region. Figure 20 
shows the densities at several altitudes for a 0.01 KT y output at 100 
km altitude, 

Having determined the secondary-electron density as a function of 
time and altitude, we can calculate the absorption of the electron 
radiation propagated to the observer, expressed in dB/km. Since the 
conductivity of a given density of secondary electrons (i.e., g/n) is 


(13) and since the 


a function of the existing electric field intensity 
field intensity is to be calculated, it is initially necessary to 

estimate the field intensity. For the estimated field value and the 
electron density values previously calculated, we can calculate the 


resulting conductivity and absorption and proceed, as will be detailed, 


to calculate the observer field. The observer's location can be either 


Fraction of secondary ionization produced, N/N, 
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Fig. 18-~ Fraction of secondary electrons from 
first Compton electron versus time 
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N,oi¢ (electrons/cc) 


-01 KT (Fig. 10) 
100-km source altitude 
Vertical geometry 


Time (shakes) 


Fig. 20— Electron density due to first-scatter Compton electrons 
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within the Compton source region or external, as illustrated in Fig. 1. 
It turns out that the field in the principal source region can be 
quite adequately approximated by simple inverse distance scaling. 

Thus, at 30 km altitude for the geometry of Fig. 1, the field will be 
about ae times the observer field. 

If we are satisfied with determination of a conservative upper- 
limit estimate of the magnitude of the observer EMP field, the numeri- 
cal complexity arising from the field-dependent conductivity problem 
can be avoided by using a constant minimum conductivity per electron, 
which is determined by the maximum field intensity value (and thus a 
minimum of absorption) at the appropriate altitude in the source region. 
After evaluation of the observer field for the assumed conductivity, 
we may refine the solution as desired by an iterative procedure, i.e., 
a new estimate of conductivity is based on the first calculated field 
and the same procedure followed with the new estimate to obtain an 
improved calculated field, and so forth. For the dB/km calculations 
of Table 2, a conductivity value 1/20th of the. value commonly used for 
low fields was assumed; hence, the dB/km values for a given electron 
density are simply 1/20th of that deduced in Ref. 14. This low value 
of effective conductivity is appropriate for electric fields of about 


(15) Hence, if the solution using this 


20,000 v/m at 20 km altitude. 
value produces indicated peak fields of, say, 5000 v/m, we have under- 
estimated absorption at 20 km altitude, and iteration of the solution 
with a properly increased conductivity will reduce the 5000 v/m 

initial calculation value. 

Inspection of Table 2, which is based on a y output of 0.01 KT, 
indicates that self-absorption is small throughout the height range of 
20 to 40 km. Thus the solution previously illustrated (Fig. 16 and 
pp. 29-31) will not be greatly affected by inclusion of self-absorption. 
For example, the signal contributed to the observer field at 4 shakes 
by electrons on the disc at a radial distance of 500 meters (2-shake 


maximum age) in Fig. 16 turns out to be attenuated less than 1 4B. 
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Table 2 


DIFFERENTIAL ABSORPTION FOLLOWING A 0.01 KT 
DETONATION AT 100 km--VERTICAL GEOMETRY 


Time After y dB per Path Kilometer 


Front Arrival 
(Shakes) 


30 40 
Altitude (km) 
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F. SUMMARY COMMENTS ON SUMMED ELECTRON-FIELD METHOD 


This section has discussed in detail the procedure for calculating 
the contribution to the observer field at any specified time from a 
thin layer in the Compton electron source region. The method of 
calculation and the thread of a numerical example were presented 
together in a step-by-step approach to make the procedure as clear as 
possible as weil as to provide a better feeling for the comparative 
aspects of the (1) delta-function versus time-variable y source, (2) 
effect of inclusions of atmospheric scatter in the calculations, and 
(3) effect of preionization, when present, on the field calculations. 
The final objective is, of course,to specify for a given geometry and 
burst parameters the electromagnetic signal amplitude versus time 
characteristics. To complete this objective, the analyst merely per- 
forms several calculations of the type given to produce Fig. 16, 
modified as appropriate to include preionization absorption, when 
present. Exactly the same procedure for finding the contribution to 
the observer field illustrated for the thin disc at 30 km is repeated 
at other altitudes in the Compton electron source region, and the 
results are integrated over height. This procedure, repeated at a 
number of observer times, thus determines the volts/meter at the ob- 
servation point as a function of time. All the numerical procedures 
required for this total solution have been brought out in this section. 
Although vertical geometry was used for ease of illustration, the 
method is identical for oblique paths if we neglect the minor 
inaccuracy brought about by the fact that the disc, as used for 
vertical geometry, will be tilted for oblique burst/observer geometry. 
Because different portions of the dise are at different altitudes, the 
circular symmetry for disc calculations used for vertical geometry is 
not strictly correct; however, for disc sizes up to a few kilometers, 
the simplifying assumption of circular symmetry about the burst/ 
observer axis is, for practical purposes, quite adequate. 

The procedure outlined in this report leads to “upper-limit" 
calculated pulse values because of omission of two factors. These are 
(1) full energy conservation considerations and (2) effects of the 


fields of the pulse itself on the orbits and hence on the radiation of 
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the individual electrons. Our work has indicated that calculations 
which omit these factors and produce calculated fields under 10,000 v/m 
are little affected (i.e., reduced in amplitude) by their inclusion; 
hence the values illustrated in this section, while upper limit in 
nature for the assumed parameters, are essentially the same as would 

be obtained by including the additional factors with their sizable 
computational complexity. Final accuracy of the result is limited 

more by basic parameters, such as conductivity, knowledge of the 

state of atmospheric ionization at the time of’ burst, and scatter 


evaluation. 
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